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ABSTRACT 

The double pulsar J0737-3039 has become an important astrophysical laboratory 
for testing fundamental physics. Here we demonstrate that the low measured 
mass of Pulsar B can be used to constrain the equation of state of neutron star 
matter under the assumption that it formed in an electron-capture supernova. 
We show that the observed orbital parameters as well as the likely evolutionary 
history of the system support such a hypothesis and discuss future refinements 
that will improve the constraints this test may provide. 
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1 INTRODUCTION 

The discovery of the first double pulsar, J0737-3039 
(Burgay et al. 2003; Lyne et al. 2004) consisting of two 
pulsars with spin periods of 22.7 ms (Pulsar A) and 2.77 s 
(Pulsar B) in a 2.4 hr orbit, has opened up a new win- 
dow for testing fundamental physics under extreme con- 
ditions: not only will the system soon allow some of the 
best tests of General Relativity (Kramer et al. 2005), 
but also it is the first system where the physics of in- 
teracting magnetospheres from two pulsars can be stud- 
ied (e.g. Arons et al. 2005). The orbit is decaying due to 
the loss of angular momentum by gravitational radiation, 
and the two neutron stars are expected to merge in only 
~ 85 Myr, a much shorter time than in any other double 
neutron star system. 

One of the other interesting characteristics of J0737- 
3039 is that the mass of the second neutron star to be 
formed (Pulsar B) is the lowest reliably measured mass 
for any neutron star to date: 1.249 ± 0.001 M0 (Kramer 
et al. 2005). Such a low mass may be an indication that 
the neutron star did not form in a standard iron-core- 
collapse supernova but in an electron-capture supernova 
(Nomoto 1984; Podsiadlowski et al. 2004) ^ These occur 

* E-mail: podsi@astro.ox.ac.uk 

^ Note, however, that it is not entirely clear that electron- 



for ONeMg white dwarfs when the core density reaches 
a critical value at which electron captures (e-captures) 
onto Mg (and subsequently Ne) start, causing a loss of 
hydrostatic support in the core and triggering its col- 
lapse. One of the key aspects of an e-capture supernova 
is that the collapse takes places when the core density 
reaches a well-defined critical value (~ 4.5 x 10^ gcm~^) 
which in turn occurs when the ONeMg core has grown to 
a well-defined critical mass (~ I.37M0; see § 3). There- 
fore, if Pulsar B indeed formed in an e-capture supernova, 
this would be the first instance for which the masses of 
the pre-collapse core and the post-collapse neutron star 
could both be determined, the former from a theoretical 
estimate of the critical mass for an e-capture supernova, 
the latter directly from the observed orbital parameters. 
Along with the pre-collapse gravitational mass, the cor- 
responding baryon number can also be calculated. Since 
the loss of material during the formation of the neutron 
star is expected to be extremely small in this scenario 
(see § 3), this is also a good approximation to the baryon 
number of the neutron star. It is the purpose of this paper 



capture supernovae necessarily produce the lowest remnant 
masses, since iron cores with masses as low as 1.15 Mq may 
be unstable to collapse (see the discussion in Woosley, Heger 
& Weaver 2002). 
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to demonstrate that both the observed orbital parame- 
ters of the system (in particular the low eccentricity) and 
the most likely evolutionary history of the system favour 
the formation of Pulsar B in an e-capture supernova and 
that comparison of its gravitational mass with the es- 
timate obtained for the baryon number enables useful 
constraints to be placed on the neutron-star equation of 
state (EoS)^. 

In § 2 we review the two most likely evolutionary 
channels that load to systems like the double pulsar and 
show how this discussion supports the key assumption 
of Pulsar B having been formed in an e-capture super- 
nova. In § 3 we provide a theoretical estimate for the 
pre-collapse core mass (with an estimate of the uncer- 
tainties) and in § 4 wo demonstrate how the properties 
of the system can be used to constrain the neutron star 
EoS. 



2 THE EVOLUTIONARY HISTORY OF 
J0737-3039 

There are two major evolutionary channels to form dou- 
ble neutron stars like the double pulsar which we shall 
refer to in this paper as the standard channel (e.g. Bhat- 
tacharya & van den Heuvcl 1991) and the double-core 
channel (e.g. Brown 1995), respectively. In the standard 
channel (see the left panel in Fig. 1), J0737-3039 origi- 
nates from a massive binary in which the more massive 
star transfers its envelope to its companion star via sta- 
ble Rochc-lobc overflow (RLOF) before it collapses in a 
supernova (SN) explosion to form the first neutron star. 
If the SN explosion does not disrupt the system, the bi- 
nary now consists of a neutron stax and a massive main- 
sequence star, and the system will evolve into a high- 
mass X-ray binary (obscrvationally it will initially look 
like a Be X-ray binary). As the secondary evolves, there 
will be a point at which it will fill its Roche lobe and 
start to transfer matter to the neutron stax. Because of 
the large mass ratio of the system, this mass transfer 
is unstable and leads to a common-envelope (CE) and 
spiral-in phase, in which the neutron star spirals towards 
the centre of the massive companion inside the compan- 
ion's envelope. Provided that the CE phase docs not load 
to the complete merger of the two components and that 
the CE is ejected, the system evolves into a very close 
binary containing the helium core of the secondary and 
the neutron star. Depending on the mass of the helium 
star and the period of the system, another phase of mass 
transfer may occur, where the neutron star is spun up 
and becomes a fast 'recycled' pulsar. Eventually the he- 
lium star collapses to form the second neutron stax in the 
system. 

In the double-core channel (sec the right panel in 
Fig. 1), the binary components are very close in mass 
initially (within 5- 10 %), and the orbit is relatively wide, 

A different test of the EoS of neutron-star matter can be 
derived from measuring the moment of inertia of Pulsar A from 
the effects of spin-orbit coupling, as proposed by Morrison et 
al. (2004). 



so that the primary only fills its Roche lobe after it has 
completed helium core burning (so-called Case C mass 
transfer) and has developed a CO core. At this stage, 
the secondary has already completed its hydrogen core 
burning phase and has evolved off the main sequence. 
Because of the high mass transfer rate, the accreting star 
expands to fill and ultimately overflow its Roche lobe, and 
the system is again expected to enter into a CE phase; 
but, in this case, the CE is formed from the combined 
envelopes of both stars. Inside the CE, there are the cores 
of the two stars, the more evolved one with a CO core 
and the less evolved He core of the secondaxy. The cores 
spiral-in inside the joint envelope until the envelope is 
ejected, leaving a very close binary consisting of the two 
evolved cores. The CO core soon collapses to form the 
first neutron star, leaving a binary consisting of a neutron 
star and a helium star. The further evolution is almost 
identical to the evolution in the standard scenario. 

Because of the constraints on the initial mass ratio 
and the orbital separation, the double-core channel is ex- 
pected to have a lower occurrence rate than the standard 
channel (by a factor of 2 to 10; Dewi, Podsiadlowski & 
Sena 2005; in preparation)^. 

The formation of J0737-3039 via the standard chan- 
nel has been studied by Dewi & van den Heuvel (2004) 
and Willems & Kalogera (2004); both studies concluded 
that this system must have originated from a close he- 
lium star-neutron star (HeS-NS) binary where the system 
underwent mass transfer during the helium-star phase, 
spinning up the first-born neutron star in the process. Be- 
cause the final stages of evolution in the standard and the 
double-core channel are essentially the same, i.e. the HeS- 
NS phase, the following discussion, which assumes the 
standard channel, also applies to the double-core chan- 
nel. 

For different reasons, Dewi & van den Heuvel (2004) 
and Willems &: Kalogera (2004) found a similar pre- 
SN mass of the helium star progenitor of J0737-3039. 
Willems & Kalogera (2004) took the threshold helium- 
star mass for the formation of a neutron star to be 2.1 Mq 
(Habcts 1986) as the lower limit. However, one should 
note that after the mass-transfer phase, the immediate 
pre-supernova mass can be as low as the Chandrasekhar 
mass (~ 1.4 M0). Dewi & van den Heuvel (2004) used 
2.3 Mq as the minimum possible helium-star mass at 
the time of the explosion, based on the assumption that 
lower-mass helium stars experience a further CE phase 
at the end of their evolution (Dewi et al. 2002; Dewi & 
Pols 2003). This lower limit on the pre-SN mass then re- 
quired a minimum kick velocity of 70 kms~^, since the 
supernova mass loss on its own would produce a much 
larger post-supernova eccentricity than is consistent with 
the present orbital parameters of J0737-3039. 

However, the assumption of the occurrence of an- 



^ However, if the spiral-in of a neutron star in a massive enve- 
lope always leads to hypercritical accretion onto the neutron 

star and its conversion into a black hole, as argued, e.g., by 
Brown (1995), the double-core channel would be the only one 
of these two channels that could produce double neutron star 
systems (also see Chevalier 1993). 
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other CE phase in a lower-mass hehum star is still an 
open question. A recent population synthesis study of the 
formation of double neutron stars (Dewi, Pols & van den 
Heuvel, 2005, in preparation) suggests that, particularly 
to explain the formation of J0737-3039, it is more likely 
that this CE phase does not occur. In this case, the final 
pre-supernova mass can be much less than 2.3 Mq, in- 
deed it can be as low as ~ 1.4 Mq, and no supernova kick 
is required to compensate for the mass loss; in particular 
it allows for the possibility that the supernova was sym- 
metric (Dewi & van den Heuvel 2005), strongly favouring 
an e-capture SN for the second supernova. 

An electron-capture supernova to form Pulseir 
B? 

The eccentricity of the double pulsar binary (c ~ 0.088 
at present, most likely 0.11-0.12 immediately after the 
second supernova; e.g. Burgay et al. 2003) is surprisingly 
low, much lower than one would expect if the system 
received a large supernova kick in the second supernova 
that formed Pulsar B. Indeed, such a low eccentricity can 
be most easily explained by a symmetric second super- 
nova in which a moderate amount of mass is expelled 
from the system. In this case the eccentricity is given by 
e = AM /{Ma + Mb), where AM is the mass lost in the 
supernova and Ma and Mb are the present masses of 
Pulsars A and B, respectively. Taking Ma = 1.338 M0 
and A^B = 1.249 (Kramer et al. 2005) and assuming a 
post-SN eccentricity eo = 0.12 then yields a prc-SN mass 
of the helium star of ~ 1.56 M©. Such low-mass pre-SN 
helium stars typically form from HeS-NS binaries with 
initial helium stars of less than ~ 3Mq (see Dewi et 
al. 2002; Ivanova et al. 2003). This includes the mass 
range where helium stars are expected to end their evo- 
lution in an e-capture supernova (Nomoto 1984). This 
is also consistent with the speculation by Podsiadlowski 
et al. (2004) that e-capture supernova© may not produce 
large supernova kicks since the explosion may proceed on 
a timcscale that is much shorter than the timescale on 
which the instabilities that produce large kicks can de- 
velop; this suggestion has received some theoretical sup- 
port from recent core-collapse calculations (Scheck et al. 
2004; H.-Th. Janka 2005 [private communication]). Since 
in any simple accretion model one would expect that the 
spin of Pulsar A would become aligned with the orbital 
momentum axis and since this orientation is not affected 
by a symmetric supernova, this may make the testable 
prediction that the post-SN misalignment angle between 
the spin of pulsar A and the orbital a^cis should be rela- 
tively small. This may indeed account for the surprising 
stability of the pulse shape of Pulsar A (Manchester et 
al. 2005). 

The second supernova kick and the space 
velocity of J0737-3039 

A low supernova kick velocity, as suggested by the low 
eccentricity, would typically also imply that the binary 
system should only have received a relatively small ad- 
ditional kick in the second supernova and that the sys- 



tem space velocity relative to the local standard of rest 
would not be much affected by it. Indeed, the space veloc- 
ity of J0737-3039 may provide an important constraint 
on its evolutionary history. Unfortunately, the situation 
regarding the value of this quantity is at present some- 
what confused. Using interstellar scintillation to measure 
the transverse space velocity of the system, Ransom et 
al. (2004) determined a large system velocity of at least 
140kms~^. Subsequently, Coles et al. (2005) showed that 
the dispersion across the field was highly variable reduc- 
ing the estimate of the scintillation velocity to a value as 
low as 66kms~^ and possibly even lower, since this value 
does not account for the motion of the Earth. More re- 
cently, Kramer et al. (2005) argued that the present limits 
on the proper motion of the system suggest a low trans- 
verse velocity of less than 30kms~'^, which would imply 
that it is statistically unlikely that the system received a 
large kick in the second supernova. On the basis of the 
original high estimate of the space velocity. Ransom et 
al. (2004) and Willems et al. (2004) concluded that the 
system should have received a large kick in the second 
supernova. However, they only considered the standard 
scenario in which the system is still fairly wide at the 
time of the first supernova; this implies that it cannot 
receive a large kick in the first supernova and remain 
bound (see Fig. 1). In contrast, in the double-core sce- 
nario where the system is already very tight at the time 
of the first supernova, the system is expected to receive 
a large systemic kick from the first supernova (of order 
150 -400 km s~^: see Fig. land Dewi et al. 2005, in prepa- 
ration) and hence no additional kick from the second su- 
pernova would be required. As far as testing the EoS is 
concerned, it is not important whether the system veloc- 
ity is low or high, since both can be understood within 
the framework of either of the two channels and a small 
second kick. In particular, we note how similar the final 
phases are in the two channels. A resolution of the issue 
of the system's space velocity could, liowcver, provide a 
powerful discriminant between the standard channel and 
the double-core channel. 



3 ELECTRON-CAPTURE SUPERNOVAE 

An electron-capture supernova occurs when the central 

density of an ONeMg core reaches the threshold value pth 
for electron captures on ^''Mg. This decreases the electron 
pressure and the electron fraction Ye , lowering the Chan- 
drasekhar mass and triggering the collapse of the core 
(Miyaji et al. 1980). 

In order to estimate the uncertainties in the criti- 
cal mass of the collapsing core, we performed a scries of 
stellar structure calculations assuming non-rotating cores 
in hydrostatic equilibrium with a prescribed central den- 
sity, homogeneous composition and a specified thermal 
profile. Since the heat released by the electron captures 
gives rise to a convective core, we adopted an isentropic 
thermal profile. 

As our reference model we adopted the composi- 
tion given by Gutierrez et al. (1996) with X(^^O)=0.72, 
X(^''Ne)=0.25, X(^^Mg)=0.03, central density pth = 
4.5 X 10^ g cm~^ and a range of central temperatures from 
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lO'^ to 10^ K. We used the equation of state of Pols ct al. 
(1995), assuming full ionization (i.e. we discard the ion- 
ization pressure term) . Note that neglecting the Coulomb 
corrections in the equation of state would increase the to- 
tal mass by 3.79X IQ-^Mq for our reference composition. 

We integrated the general relativistic (GR) equations 
of hydrostatic equilibrium out from the centre 




Here m and A are the gravitational mass and baryon 
number enclosed within a sphere of radius r; p is the 
density, including contributions from both the rest mass 
and the thermal energy; Ub is the baryon number den- 
sity. Rather than talking in terms of the baryon num- 
ber A, which is a rather abstract quantity, it is con- 
venient to convert this into a mass by multiplying by 
the atomic mass unit (931.50 Mev/c^), and we refer to 
this quantity as the haryonic mass. If we had ignored the 
general relativistic corrections, our estimate of the bary- 
onic mass would have boon increased by 1.30 x IO^^Mq 
(~ 1%). This difference is larger than might have been ex- 
pected; the reason, however, is analogous to that for the 
large difference between the Newtonian and GR values 
for the maximum mass of a neutron star (Oppenheimer 
& Volkoff 1939). 

To check the validity of our procedure, we com- 
puted a model without the GR corrections to repro- 
duce the case considered by Miyaji et al. (1987) and 
Nomoto (1987), where they used a composition X(^^0) = 
0.12, X(2°Nc) = 0.76, xl^%lg) = 0.12, central density 
3.98 X 10^ gcm~^ and central temperature log(rc) = 8.61. 
We find a total mass m = 1.3754 M©, which is very close 
to their published value of 1.375 M©. 

This composition is, however, no longer considered 
appropriate since the reaction rate ^^C(a, 7)^^0 has been 
revised upwards (Dominguez, Tornambe & Isern 1993). 
This leads to a lower Ne and Mg abundance and a higher 
O abundance at the end of C burning than computed 
with the Fowler et al. (1975) rates (as was done in Miyaji 
et al. 1987). To estimate the uncertainties introduced by 
the composition, we investigated several different com- 
positions found in the literature (Dominguez et al. 1993; 
Gutierrez ot al. 1996; Ritossa, Garcia-Borro & Ibon 1996; 
Gil-Pons & Garcia-Bcrro 2001) after the revision of the 
^^Claj^y^'O rate. The composition is important in de- 
termining the Coulomb parameter and Ye, and hence the 
magnitude of the Coulomb corrections. We use our ref- 
erence model to give a lower bound on the critical mass 
(solid curve in Fig 2) and the composition X(^^0) = 0.56, 
XC°Ne) = 0.29, X(2*Mg) = 0.06, ^(^^Na) = 0.07 and 
X{^^C) = 0.01, which mimics the composition of Gil- 
Pons & Garci'a-Berro (2001) and has the lowest Coulomb 
correction, to give an upper bound to the critical mass 
(dashed curve in Fig 2). 



Finally, wo changed the threshold density pth from 
4.5 X lO'^ to 4 X 10^ g cm ^ in order to estimate the effect 
of the shift in the critical density due to the Coulomb 
corrections (the latter is the appropriate density without 
any Coulomb corrections; see Gutierrez ct al. 1996 for 
a detailed discussion). This change decreases the critical 
mass by 2 X 1O~^M0 (dotted and dot-dashed lines in 
Fig 2). 

1) In the previous studies, the central temperature after 
the onset of electron captures ranged from log Tc = 8.5 
up to log Tc — 8.65. For this temperature range. Figure 2 
then implies that the baryonic mass of the pre-coUapse 
core should lie between 1.366 and 1.375 Mq (the thin 
dotted curves) taking into account the uncertainties in 
temperature and composition. On the assumption that 
the loss of material during formation of the neutron star 
is negligible, this then gives the predicted range for the 
baryonic mass of the neutron star (which we refer to as 
Mo). 

Caveats 

At this point, several caveats should be made about 
other effects which may systematically affect this esti- 
mate. First, we take the baryonic mass of the pre-collapse 
ONeMg core to be the same as that of the neutron star, 
neglecting any loss of material during the formation of the 
neutron star. In practice, some material may be ejected in 
the supernova in a neutrino-driven wind (Qian & Woosley 
1996). However, because of the steep density gradient at 
the edge of the ONeMg core, we expect this mass loss 
to be small, probably less than a few times 10""^ M©, al- 
though it could potentially be as large as 10~^ Mq and 
the amount of mass loss itself depends on the EoS (Janka 
2005, private communication). Second, an e-capture su- 
pernova may occur between carbon shell flashes when the 
central density increases significantly (although we note 
that this point is not yet fully resolved). Because of the 
discrete nature of the carbon fiashes, this may introduce 
a natural variation in the critical mass from star to star 
because of variations of the thermal profile in the outer 
ONeMg core. In this context we note that the assump- 
tion of an isentropic profile is unlikely to be correct in 
the outer parts of the ONeMg core. 

It is clear that these effects could significantly change 
our estimate for the baryonic mass of the neutron star 
(probably decreasing it) and therefore our present esti- 
mates should only be considered as preliminary. How- 
ever, with the expected progress in simulating e-capturc 
supernovae and calculating the evolution of their pro- 
genitors (with improved e-capture rates, a richer nuclear 
network, inclusion of accurate Coulomb corrections and 
using GR), we estimate that ultimately one should be 
able to pin-point the mass of the collapsing core to within 
~ 2 X IQ-^Mq. 



4 CONSTRAINTS ON THE EQUATION OF 
STATE 

On the hypothesis that the scenario presented in the 
previous sections is correct, we can use the information 
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Figure 2. Baryonic mass versus central temperature for vari- 
ous ONeMg cores in hydrostatic equilibrium. Solid curve: ref- 
erence model with p^jj = 4.5 x 10^ gcm~^ with O and Ne mass 
fractions X(^^0) = 0.72 and X(^°'Ne) = 0.25 (minimum mass 
for the variation in composition). Dashed curve: similar but 
with X(l'50) = 0.56 and X(^"Nc) = 0.29 (maximum mass for 
the variation in composition). Dotted curve: pth decreased to 
4 X 10''gcm^'^ with the composition of the reference model. 
Dot-dashed curve: same model as for the dashed curve but 
with decreased pth- The thin dotted lines give the most likely 
range of parameter values at the time of core collapse. 



about the gravitational and baryonic masses of Pulsar B 
to place constraints on the EoS of neutron-star matter. 
These turn out to be quite interesting. 

For any given EoS for neutron-star matter, one can 
calculate the relation between the gravitational mass and 
the baryonic mass (bearing in mind that the rotation 
speed is so low that taking the object to be spherical and 
non-rotating is an excellent approximation). The present 
gravitational mass (known from observations) and the 
baryonic mass (known from the stellar evolution calcular- 
tions) then define an error box through which the rela- 
tions calculated from the equations of state would need to 
pass in order to be consistent. In this section, we discuss 
how the constraints obtained in this way turn out. The 
observed gravitational mass {Mq = 1.249 ± 0.001 M©) 
and the calculated baryonic mass (Mo in the range 
1.366—1.375 Mq) specify the boundaries of the error box. 

Results obtained by integrating the GR equations of 
hydrostatic equilibrium (Eqs 1-3) for a range of EoSs are 
shown in the four panels of Figure 3, each corresponding 
to a particular class of equations of state (we use a mod- 
ified form of the classification in the paper by Morrison 
et al. 2004) . All of these EoSs give the maximum gravita- 
tional mass Afnmx for a non-rotating neutron star as being 
above 1.5 Mq, in line with observations. Class I EoSs (top 
left-hand panel) come from non-relativistic many-body 
calculations with "realistic" potentials: APR98 (Akmal, 
Pandharipande & Ravenhall 1998), WFF88 (Wiringa, 
Fiks & Fabrocini 1988) and EPS (Lorenz, Ravenhall & 
Pethick 1993) include only nucleonic degrees of freedom, 
BJ74 (Bcthc & Johnson 1974) and MOSZ74 (Mozskowski 
1974) include also hyperonic components at the higher 



densities. Class 11 EoSs (top right-hand panel) use rola- 

- tivistic mean-field (or effective-field) approximations in- 
cluding hyperonic degrees of freedom: GLE210, GLE240 
and GLE300 (Glendenning 2000), and HOFOl (Hofmann 
et al. 2001) arc shown. We also include hero one EoS 

''Z representing a hybrid stellar model (nucleons -I- quarks): 
GLENHYB (Glendenning 2000). Class III EoSs (bot- 
tom left-hand panel) use non-relativistic phenomenolog- 

- ical potentials of the Skyrmc type (see Stone et al. 2003 

- and references therein). Class IV EoSs (lower right-hand 

- panel) are for other phenoincnological non-relativistic po- 

- tentials: BPAL21 and BPAL31 (Prakash et al. 1997) have 

- only nucleonic degrees of freedom, while BAL97 (Bal- 

- berg & Gal 1997) includes hyperons at high density. All 
: of these high-density EoSs were joined onto the Baym- 

^jgBethe-Pethick EoS (Baym, Bethe & Pethick 1971) at a 
density of ~ 1.4 x lO" g cm"^ (0.08 - 0.09 ftn"^) and, in 
turn, this was joined onto the Baym-Pethick-Sutherland 
EoS (Baym, Pethick & Sutherland 1971) at 4.2 x lO" 
gcm~^ (2.5 X 10~*fm~^). By doing this, the inner and 
outer crust of the neutron star were treated in the same 
way for all of the EoSs, and so all of the differences 
seen result from differences in the treatment of the high- 
density matter. (Strange star models have not been in- 
cluded in Fig 3; they cover a very wide range and give 
curves passing both above and below the error box as 
well as curves passing through it.) 



The most clear-cut result is that none of the class II 
models tested in this work give predictions in line with 
our constraint. For the other classes (I, III and IV), the 
situation is less clear and depends on the particular prop- 
erties of each individual EoS. For the phenomenological 
Skyrme potentials, the EoSs give a wide range of pre- 
dictions within the region delimited by Skll and BSk8. 
Those parametrisations giving Mq/Mq curves passing 
through the error box give Mmax between 1.6 — 1.9 Mq. 
All of those for which Mmax > 1.9 Mq give curves passing 
above the error box while those for which Mmax < 1-6 Mq 
give curves passing below it. With reference to the dis- 
cussion in Stone et al. (2003) and using the notation of 
that paper, we note that all of the Skyrme EoSs pass- 
ing through the error box are type II parametrisations 
whereas all of those coming from type I parametrisations 
pass above it. 



Apart from the comments made above, there is no 
simple general interpretation of the implications of our 
constraint for the physics behind the particular EoSs. It 
is important to recognise that our constraint (assuming 
that our overall scenario is correct) represents a neces- 
sary but not sufficient condition for choosing a suitable 
EoS for neutron star models. Additional observational 
information is needed, in particular concerning the neu- 
tron star maximum mass which, in combination with the 
present constraint, would give a more definitive criterion 
for the choice of physical model for the EoS. Also, the 
influence of different treatments for the inner and outer 
crust needs to be investigated before any final conclusion 
is drawn. 
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Figure 3. Relation between the gravitational mass Mq of neutron star models and their baryonic mass Mq, measured in units 
of the solar mass Mq , for various equations of state. The constraint derived in this paper is marked by a rectangle. For more 
explanation, see text. 



5 CONCLUSIONS 

In this paper, we have demonstrated that the measured 
gravitational mass of Pulsar B in the double pulsar 
J0737-3039 can be used to give a new test of the neutron- 
star EoS if one makes the all-important assumption that 
Pulsar B was formed in an electron-capture supernova. 
In this case, its baryonic mass can be estimated theo- 
retically and comparison between this and its gravita- 
tional mass can then be used to constrain the EoS. We 
have re-constructed the possible evolutionary histories 
for J0737-3039 in the main formation channels to sup- 
port the hypothesis of Pulsar B having originated in an 
electron-capture supernova and have discussed possible 
tests of this hypothesis. Future refinements, both of the 
stellar evolution models (to better pin down the critical 
pre-coUapse mass) and of electron-capture collapse mod- 
els (to quantify the possible mass loss), should lead to 
a more stringent constraint which, when combined with 
other astrophysical EoS constraints, can provide new in- 
sight into the physics of neutron-star matter. 
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